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Abstract: Quasi-simultaneous optical/near-IR photometry is presented for a 
sample of 37 luminous optically selected QSOs drawn from the Large Bright 
G> ' QSO Survey. Most of the QSOs have decreased in brightness since discovery; 

^ ' this is expected in flux-limited samples. The continuum shape of most of the 
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QSOs can be represented by a power-law of the form F(iy) oc u , but a few 



5^ ■ have softer (redder) continuum slopes. 
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1. Introduction 

The shape of the continuum spectra of AGN in the rest-frame UV-optical 



region has long been uncertain. It is typically parameterised as a power-law 
of index a, where the flux per unit frequency f^, oc 1'°'. The mean value 
of a, the dispersion in values of a, and even the validity of this power-law 
^ parameterisation, are all controversial. 

■ For many years, a mean value of a = —0.7 has been assumed by many 

researchers, based largely on the observations of Richstone & Schmidt (1980) 
(see also O'Brien, Gondhalekar & Wilson 1988 and Sargent, Steidel & Bok- 
senberg 1989). Recently, however, a number of authors have claimed that the 
true mean continuum slope is actually much harder (bluer) than this, with 
a = —0.3 (eg. Neugebauer et al. 1987, Francis et al. 1991). There is equal 
disagreement about the range of continuum slopes shown by QSOs, with some 
authors claiming that all QSOs have essentially the same continuum slope (ie. 
little or no dispersion in a, Sanders et al. 1989) while others claim that the 
dispersion in continuum slopes is large (Elvis et al. 1994, Webster et al. 1995). 

These disagreements are important both physically, and because of the 
enormous uncertainty they introduce into measurements of QSO evolution. 
If, for example, the mean value of a is —0.3, rather than —0.5 as commonly 
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assumed, the number counts of bright QSOs ai z ^ 2 can be explained with 
only half the evolution in the comoving density of luminous QSOs previously 
required, due to the difference in /c-corrections. Allowing for a reasonable 
scatter in continuum slopes can also introduce factors of two or more into 
the inferred evolution of QSO number densities (Giallongo & Vagnetti 1992, 
Francis 1993b). 

Physically, a value of a ~ —0.3 is consistent with the free-free emission 
models of Barvainis (1993), whereas a large dispersion in continuum slopes, 
unless caused by dust, would rule out this model. Webster et al. (1995) 
claim that fiat-spectrum radio-loud quasars show a much greater dispersion in 
their continuum slopes, and a much softer (redder) mean slope than optically 
selected QSOs, and ascribe this difference to dust obscuration, a claim disputed 
by Serjeant & Rawlings 1996, who ascribe the results to relativistically beamed 
red synchrotron emission. 

The different values of a reported in the literature may reflect real differ- 
ences between different classes of QSOs, such as radio-loud and radio-quiet, or 
luminous and faint. It is however possible that the differences reflect no more 
than different measurement techniques. One problem is that most AGNs with 
multi-wavelength observations are low luminosity, so the contribution of the 
host galaxy light is important (Elvis et al. 1994) and needs to be carefully 
corrected for. Measurement in the optical is complicated by the many strong 
blended Fe II emission-hnes that obscure large parts of the spectra (Wills, 
Netzer Sz Wills 1985, Francis et al. 1991). This obscuration makes measure- 
ments of continuum slopes inferred from optical data alone very unreliable, a 
problem that can be solved by the use of multi-waveband data. Alas, such 
data is seldom taken simultaneously, allowing variability to affect the inferred 
continuum slopes. 

In this paper, an attempt is made to provide a definitive measurement of 
the continuum slope distribution of one type of quasars; luminous optically 
selected QSOs. The QSOs are drawn from the Large Bright QSO Survey 
(LBQS, Hewett, Foltz & Chaffee 1995), a large, uniform and highly complete 
optically selected QSO sample that is increasingly being used as a standard 
against which to compare other quasar samples (eg. Webster et al. 1995). 
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These QSOs are bright enough that the emission from their host galaxies can 
be neglected, as shown below. To minimise the effect of blended emission-lines, 
the continuum slopes are inferred from multi-waveband optical and near-IR 
photometry covering a factor of five in wavelength. The photometry is quasi- 
simultaneous, to minimise any variability effects. This data set will provide 
a benchmark against which the continuum slopes of other samples can be 
compared. 

The sample and observations are described in Section 2, the results are 
presented and discussed in Section 3, and conclusions drawn in Section 4. 

2. Observations 

2a. The Sample 

A subset of QSOs drawn from the Large Bright QSO Survey (LBQS) were 
observed. The LBQS is described in Foltz et al. (1987, 1989), Hewett et al. 

(1991) , Chaffee et al. (1991) and Morris et al. (1991), and summarised in 
Hewett, Foltz and Chaffee (1995). It consists of 1065 QSOs, with redshifts 
0.2 < z < 3.4, apparent magnitudes 16.0 < Bj < 18.85, and absolute mag- 
nitudes Mbj < -21.5 (for go = 0.5, Hq = 100km s'^Mpc"^). The QSOs 
were selected from UK Schmidt objective prism plates, using a variety of au- 
tomated selection criteria, including continuum colour and the presence of 
strong emission-lines and continuum breaks. Comparison with existing sur- 
veys suggests that the LBQS is highly complete to its Bj magnitude limit. 
Radio observations of a sub-set of the LBQS are presented by Visnovsky et al. 

(1992) , Francis, Hooper & Impey (1993) and Hooper et al. (1995). Discovery 
epoch Bj magnitudes were obtained by scans of UK Schmidt plates and are 
accurate to ~ 0.15 mag. 

B, V, H and K photometry for subset of 37 LBQS QSOs is presented 
in this paper. The subset observed were selected to be optically bright, high 
redshift, and accessible during the observing runs. 

At redshifts z ~ 2, typical of most of the subsample observed, the iC-band 
corresponds to the rest-frame it! band, and the S-band corresponds to rest- 
frame 1500A. We are therefore measuring the slope of the long- wavelength end 
of the 'Big Blue Bump' (Sanders et al. 1989). 
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The spectral synthesis models of Bruzual & Chariot (1993), and of Rocca- 
Volmerange & Guiderdoni (1988), were used to estimate the colours of the 
most luminous plausible host galaxies; starbursting giant proto-ellipticals. The 
rest-frame rest-frame i?-band magnitudes at 2 ~ 2 were be calculated, and 
converted into observed K-hand magnitudes. Even the most extreme galaxies 
had K > 17.5; two magnitudes fainter than the QSOs in this sub-sample 
(Table 1). Host galaxy light is thus unlikely to contribute significantly to 
these observations. 

2b. Observations 

Optical photometry, in the B and V bands, was obtained on the Steward 
Observatory 90" telescope on the night of 3rd September 1992. Conditions 
were photometric, but due to malfunctions in both the autoguider and the 
tracking, fast manual guiding was required. Five minute CCD exposures were 
obtained. Not all objects could be imaged, due to substantial overheads asso- 
ciated with the manual guiding. 

Near-IR photometry, in the H and Ks bands, was obtained on the nights 
of 27th and 28th April 1992, and 27th - 30th September 1992, on the Steward 
Observatory 61" telescope, using the 256 x 256 NICMOS HgCdTe array cam- 
era. 3x3 raster patterns were used, with three 50 second exposures at each 
grid pointing in H, and five 30 second exposures in Kg to give an exposure 
time of 22.5 minutes. Conditions were photometric throughout. 

The data was reduced using standard IRAFQ procedures. Magnitudes were 
measured using very large fixed apertures, as sky noise was not the dominant 
uncertainty, and the point spread function was very irregular due to the manual 
autoguiding. 

Only the September data was used to derive continuum slopes; the optical 
and IR photometry are separated by three weeks, and so are not strictly simul- 
taneous. Optically selected QSOs, however, show little variability of timescales 
of less than a month (Cristiani et al. 1996). Magnitude errors in both the op- 
tical and near-IR are ~ 0.05 magnitudes. 

^IRAF is distributed by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc. (AURA) under 
cooperative agreement with the National Science Foundation. 
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3. Results and Discussion 

The results of the photometry are shown in Table 1 . 
3a. Variability 

The UK Schmidt sky survey plates from which the original magnitudes 
were measured were taken in the late 1970s and early 1980s. The time interval 
between the original Bj magnitudes and our B and V photometry is thus ~ 10 
years, allowing the assessment of the optical variability of the sample. 

1992 epoch Bj magnitudes were calculated using the approximate relation- 
ship Bj = 0.7i? + 0.3y, and the change in Bj between the Schmidt plate epoch 
and 1992 calculated. A histogram of the changes is shown in Figure 1. 

As Figure 1 shows, the magnitude changes arc very lopsided; of the 26 
QSOs measured, 22 decreased in brightness over the ~ 10 years, while only 4 
increased in brightness. The situation is even more lopsided when larger mag- 
nitude changes are considered; of the 16 QSOs which varied by more that 0.1 
magnitudes, none increased in brightness. The median change in brightness is 
a decrease of 0.17 magnitudes, and the standard deviation is 0.23 magnitudes. 
This standard deviation is in good agreement with other studies (eg. Cristiani 
et al. 1996, and references therein). This lopsidedness remains regardless of 
the exact photometric techniques employed. 

Why are the magnitude changes so lopsided? If QSO variation is a random 
walk process, the probability of a QSO increasing and decreasing in brightness 
will be equal. If this were true, the probability of getting such a lopsided 
distribution by chance is 0.6%. 

If, however, the variability of a QSO is not a random walk process, but each 
QSO has some characteristic brightness around which it varies, the lopsided- 
ness can be explained. Any sample will consist of QSOs that were abnormally 
bright at the time of selection, and those that were abnormally faint. The 
QSO luminosity function is, however, very steep, so a QSO detected at a given 
magnitude is far more likely to be a fainter QSO caught in a bright state, than 
a brighter QSO caught in a faint state. 

The size of the expected lopsidedness can be crudely estimated as follows. 
The typical variation is around 0.2 mag. Given the luminosity function of 
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bright QSOs (eg. Boyle, Shanks and Peterson 1988), QSOs 0.2 magnitudes 
brighter than a given magnitude are ~ x4 less common than those 0.2 mag- 
nitudes fainter than the given magnitude. Thus the number of QSOs at the 
given magnitude which were at a maximum at the selection epoch is likely 
to exceed that which were at a minimum by a factor of ~ 4. Thus the ratio 
of the number QSOs that subsequently grow fainter to the number that grow 
brighter should be ~ 4:1, in rough agreement to that observed. 

3b. Continuum Slopes 

Continuum slopes were computed as follows. Only QSOs with both optical 
and near-IR photometry were used. The optical magnitudes were corrected 
for the presence of strong broad emission-hnes in the B and V bandpasses. 
For each spectrum, the redshift determined which strong emission-lines would 
be in each filter bandpass. If a strong line was present, its equivalent width 
in that QSO was obtained from the measurements of Francis (1993a). The 
flux in the filter bandpass was then corrected for the presence of this hne, 
allowing for the redshift dependence of observed-frame equivalent widths. For 
the near-IR bandpasses, which often contained H-a and H-/?, no measurements 
of the equivalent widths of individual QSOs were available, so the fluxes were 
corrected by assuming a mean rest-frame equivalent width of SSA for H-/3 and 
290A for H-a (Francis et al. 1991). 

For the higher redshift QSOs, Ly-a forest absorption can reduce the B- 
band flux. At z > 3, the whole S-band is affected by Ly-a absorption, and 
no correction was attempted; the 5-band flux was simply not used in the flt. 
For 2A < z < 3 however, the S-band flux was corrected, the correction factor 
being chosen from visual inspection of the spectrum. Magnitude corrections 
(corrected magnitude minus observed magnitude) are shown in Table 2, and 
the corrected continuum fluxes in Table 3. Note that these corrections are 
typically small; even if no corrections are employed, the results of this paper 
are almost identical. 

No attempt has been made to correct for weak emission-lines, and for Fe II 
emission. The use of a mean equivalent width for the Balmer line correction 
also introduces an error of ~ 30% in the H and K-band correction factors 
(Francis 1993a). However, the long wavelength baseline (0.4 - 2.2//m) means 
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that these errors will not have a large effect on the derived continuum slope a; 
a 25% error in one of the continuum fluxes introduces only an error of ~ 0.1 in 
the continuum slope index a. Errors in a are predominantly systematic and 
hard to estimate, but from the self consistency of objects for which we have 
multiple observations, and the residuals to power-law fits to the continuum 
(Table 3), we estimate them to be of order 0.1. 

A power-law continuum is fitted (by least-squares minimisation) to the 
corrected continuum fluxes; the resultant power-law indicies a are listed in 
Table 3. A histogram of these continuum slopes in shown in Figure 2. The 
power-law flt is generally a good one; the scatter around the flts (Table 3) is 
consistent with the photometry and correction errors. Note however that in 
the absence of R or /-band photometry, we have little sensitivity to large-scale 
departures from a power-law continuum shape. 

The mean value of the continuum slope index a is —0.46, with a root- 
mean-squared (rms) dispersion of 0.30. Thus the mean slope is harder than 
the canonical a = —0.7, but not as hard as the a = —0.3 claimed by Francis 
(1993b). Note however that the median slope is —0.35 and the mode ~ 0.3; 
the mean being biassed towards higher values by the two outlying points. 
There are no obvious observational problems with the observations of these two 
outlying points; their redness is much too large to be explained by errors in the 
photometry and/or correction factors. Interestingly, both the outlying points 
are radio-loud QSOs, though the sample is too small to draw any significant 
conclusions. 

These results demonstrate that the continuum slopes do show a real dis- 
persion, but a much smaller one than that claimed by Webster et al. (1995) 
and Elvis et al. 1994). A standard deviation of 0.3 in a agrees with that 
found by Giallongo & Vagnetti (1992), who showed that it leads to important 
consequences for measurements of QSO luminosity functions. The hard (blue) 
mean continuum slope reported here also suggests that previous workers, us- 
ing a — —0.7, will have overestimated the evolution of the number density 
of QSOs. A dispersion in a of 0.3 is smaller than that reported by Francis 
(1993b) for the LBQS as a whole (0.5); the difference is probably that the 
measurements here are more accurate, being based on a much wider wave- 
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length coverage. 
4. Conclusions 

The very hard (blue) modal continuum slope of this sample agrees well with 
the measurements of Sanders et al. (1989), for another sample of luminous, 
optically selected QSOs. It is however, very different from the continuum slopes 
of the radio-selected sample of Webster et al. (1995) and the heterogeneous 
sample of Elvis et al. (1994), both of which show a substantial dispersion in 
continuum slopes, with —0.3 > a> —4. It is interesting to note, however, that 
the continuum slopes of the optically selected QSOs lie at the blue envelope 
of the continuum slope distribution of Webster et al.'s radio selected quasars. 

Webster et al. suggested that the redness of their data was caused by dust 
obscuration, somewhere along the line of sight, and that the true continuum 
shape of QSOs is a uniform f{i/) oc v~^-^. They showed that if this were 
true for radio-quiet QSOs, optically magnitude limited samples would show a 
distribution of continuum slopes with a narrow peak at a = —0.3, and a low 
tail to redder slopes. This is a good description of our data (though given the 
small sample size, many equivalently good descriptions are possible). 

Our observations are therefore consistent with a picture in which quasars 
have an intrinsic continuum slope of a ~ —0.3 (as predicted by the free-free 
emission model), but are reddened by various amounts of dust, as suggested 
by Webster et al. Most of the heavily reddened quasars will be missed from 
samples, such as the LBQS, with a blue magnitude limit, but a few moderately 
red sources will remain, as observed. This model will be further elaborated in 
Francis et al. (1996, in preparation). 
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Table 1. Observed Magnitudes 



Name B V Old H Redshift Radio 

Bj Loud? 



0004+0147 


18 


277 


17 


672 


18 


13 


15 


99 


15 


56 


1 


7110 


? 


0006+0200 


17 


968 


17 


605 


16 


51 


16 


35 


15 


62 


2 


3358 


No 


0006+0230 


18 


569 


18 


339 


17 


99 


16 


83 


16 


00 


2 


0981 


No 


0007+0142 


18 


344 


18 


154 


18 


26 


16 


61 


16 


08 


1 


7654 


Yes 


0009+0219 


18 


360 


17 


960 


17 


99 


15 


91 






2 


6416 


No 


0010-0012 


19 


036 


18 


713 


18 


46 










2 


1536 


? 


0013-0029 


18 


436 


18 


048 


18 


18 










2 


0835 


No 


0018-0220 


17 


486 


17 


105 


17 


44 






15 


46 


2 


5960 


No 


0049+0045 


17 


766 


17 


540 


17 


50 


15 


74 






2 


2647 


No 


0052-0058 


18 


347 


18 


194 


17 


93 


16 


29 


15 


35 


2 


2120 


No 


0056-0009 


17 


929 


17 


664 


17 


73 


15 


90 









7175 


Yes 


0106+0119 


18 


498 


18 


171 


18 


34 


15 


12 


14 


04 


2 


0989 


Yes 


0256-0000 


18 


876 


17 


606 


18 


22 


15 


72 


15 


19 


3 


3638 


Yes 


0256-0031 










17 


59 


15 


49 


14 


80 


1 


9951 


? 


0257-0254 










17 


60 


14 


71 


14 


51 


1 


0700 


? 


0258+0210 


18 


120 


17 


907 


17 


97 


15 


88 


15 


52 


2 


5239 


No 


0301+0015 










18 


31 


16 


44 






1 


6458 


? 


0302-0019 


18 


409 


17 


436 


17 


78 


15 


87 


15 


27 


3 


2814 


No 


1012+0213 










17 


63 






15 


10 


1 


3778 


Yes 


1013+0124 










16 


52 






14 


48 





7790 


Yes 


1215+1121 










16 


61 






14 


96 


1 


3984 


Yes 


1216+1754 










18 


07 






15 


23 


1 


8097 


Yes 


1222+1433 










17 


11 






14 


85 


1 


3368 


No 


1233+1524 










18 


06 






16 


20 


1 


5441 


No 


2230+0232 


18 


387 


18 


038 


18 


05 


16 


43 


15 


72 


2 


1468 


No 


2231-0015 


18 


096 


17 


394 


17 


53 


15 


46 


15 


02 


3 


0150 


No 


2231-0048 


17 


623 


17 


241 


17 


57 










1 


2095 


? 


2231-0212 


18 


357 


18 


065 


18 


15 






15 


80 


1 


9052 


? 


2239+0007 


18 


796 


18 


453 


18 


28 


16 


32 






1 


4401 


? 


2241+0014 










17 


57 


17 


72 


15 


53 


2 


1310 


No 


2241+0016 


18 


361 


17 


996 


18 


30 


16 


52 


15 


74 


1 


3936 


? 


2243+0141 


18 


911 


18 


544 


18 


25 


16 


52 






2 


3140 


No 


2244-0105 


18 


091 


17 


680 


17 


95 


15 


91 


15 


29 


2 


0300 


No 


2244-0208 


19 


421 


19 


016 


15 


.4 


17 


55 






1 


9682 


? 


2247+0135 


17 


945 


17 


378 


17 


77 










1 


1284 


? 


2248+0127 


18 


582 


18 


163 


18 


21 






15 


81 


2 


5585 


No 


2350-0045B 










18 


38 






15 


69 





4437 


? 
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Table 2. Magnitude Correction Factors 



Name B V H K 



Correction to Observed Magnitude 



0004+0147 


0, 


.06 


0.03 


0.0 


0.0 


0006+0200 


0, 


.10 


0.10 


0.08 


0.28 


0006+0230 


0, 


.13 


0.05 


0.08 


0.27 


0007+0142 


0, 


.14 


0.03 


0.0 


0.0 


0009+0219 


0, 


.08 


0.11 


0.0 




0018—0220 


0, 


.12 


0.07 




0.0 


0049+0045 


0, 


.07 


0.09 


0.08 




0052-0058 


0. 


.30 


0.0 


0.08 


0.27 


0056-0009 


0, 


.08 


0.0 


0.0 




0106+0119 


0, 


.19 


0.0 


0.08 


0.27 


OZOD — 0000 






O.o 1 


0.0 


ci no 
O.Oo 


0258+0210 


0, 


.04 


0.13 


0.09 


0.0 


0302-0019 






0.48 


0.0 


0.08 


2230+0232 


0, 


.13 


0.0 


0.08 


0.27 


2231-0015 






0.0 


0.0 


0.08 


2231-0212 


0, 


.06 


0.04 




0.08 


2239+0007 


0, 


.06 


0.0 


0.28 




2241+0016 


0, 


.06 


0.0 


0.28 


0.0 


2243+0141 


0, 


.18 


0.10 


0.08 




2244-0105 


0, 


.04 


0.06 


0.0 


0.26 


2244-0208 


0, 


.09 


0.11 


0.0 




2248+0127 


0, 


.20 


0.12 




0.0 
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Table 3. Corrected continuum fluxes 



IN cxLlLxd 


u.44yLtlIl 


0.55/xm 1.6;um 
(flux X10-29W m 


2. l/iin 
n.z 


rvt:blCIU.cilb l^lllib j 

TTl 1 1 


Slope 
cx 




n on^ 
u.zuo 


n oo/i 
U.zy4 


n A QQ 


n 41 1 

U.4i 1 


n n'^f; 

U.UOO 


— U.O / / 


nnnfij-nonn 
uuuD-ruzuu 
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Figure Captions: 

Figure 1 — The distribution of magnitude changes of QSOs, between the old, 
UK Schmidt Bj magnitudes and those presented in this paper. 

Figure 2 — The distribution of continuum slopes, as measured by the power- 
law index a, where the continuum flux per unit wavelength f{u) if fit by a 
function of the form /(i/) oc u"'. 
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